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ABSTRACT 


A fheoreHcai investigation was performed to develop a model for the radar 
return from bare ground. The validity of the model was tested by comparing its 
theoretical prediction with measured data collected by the University of Kansas 
Remote Sensing Laboratory S-IS GHz radar spectrometer system. 

It was assumed that the target area consists of a collection of small, medium 
and large size facets. Then this model was used to calculate the radar cross section 
of bare ground and the effect of the frequency averaging on the reduction of the 
variance of the return . 

It was shown that by assuming that the distribution of the slope to be Gaussian, 
and by assuming that the distribution of the length of the facet to be in the form of the 
positive side of a Gaussian distribution, the results are in better agreement with 
experimental data than the results of previous facet models. It was also shown that 
for this calculation we do not need to know the exact correlation length of the small 
structure on the ground, instead an effective correlation length was calculated based 
on the facet model and the wavelength of the incident wave. Hence, the parameters 
necessary to specify the surface are: standard derivations of slope in x and y directions, 
standard deviation of the distribution of the facet size, and the dielectric constants 
of the target. 

For investigating the effect of the frequency averaging we expanded the 
previously available results based on the uniform scatterer model and took into 
consideration the penetration effect. It was shown that at small incidence angles, 
the number of independent samples predicted is significantly larger and in better 
agreement with measured data from alfalfa. 

It was also shown that based on the facet model assumption the reduction in 
the variance of the return is not only a function of the product of the sweep band 
and the time span of the target, as the uniform scatterer model indicates, but it is 
also a function of the geometrical properties of the surface, center frequency of 
the incident wave and the polarization. 
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1.0 INTRODUCTION 


The purpose of i"his invesUgafion is to develop a theoreHcal model for fhe 
radar refurn from bare ground. The validity of the model is to be tested by comparing 
its theoretical predictions with measured data collected by the University of Kansas 
Remote Sensing Laboratory 8-18 GHz radar spectrometer system. 

In 1956 Kertzin [7] proposed that the surface of an area-extended target can be 
assumed to be composed of a collection of independent facets of varying sizes and 
slopes. This facet model was developed for predicting sea backscatter at large 
angles of incidence. Over its range of validity (>85°), Katzin reports good agreement 
with experimental data. Using the same facet backscatter approach, Katzin later 
developed a second model [8] to cover the other end of the incidence angle range 
(0° to 40°). In both models Katzin divides his facet sizes into small and large facets 
and assumes certain facet size and facet slope distributions. 

A study of the backscatter from circular disks Indicates that the variation of the 
radar cross section of a disk as a function of its radius (measured in wavelengths) suggests 
that facets should be separated into three different size ranges rather than two. Thus, 
in the present study a medium size facet group is introduced along with other modifications 
of Katzin's models. The results indicate a significant improvement over Katzin's two 
scale model in terms of agreement with experimental data. Furthermore , the proposed 
model does not have any inherent incidence angle range limitations. 

Prior to developing the model proposed in this study (Chapter 3), a detailed 

discussion of Katzin’s models is presented in Chapter 2 in order to facilitate the 

understanding of the modifications made in Chapter 3. Chapter 4 describes the 

properties of the 8-18 GHz radar spectrometer used to collect the experimental data. 

Two types of experiments were performed; the first experiment consisted of making about 

140 independent backscatter measurements of different illuminated cells of alfalfa and 

bare~soil fields. Measurements were made at 9 GHz with W polarization at incidence 

angles of 0° (nadir), 10°, 20°, 30°, 50° and 70°. At each angle, the mean and 

variance of the backscattering coefficient were calculated and used to calculate the 

2 

number of independent samples assuming an X fading distribution. In the second 
experiment, the radar return from a bare ground field was measured under several soil 
moisture conditions. Along with the moisture changes, the surface texture changed 
as well . Comparison of theoretical predictions with measured data is presented in 
Chapter 5. 



2.0 TECHNICAL REIVE W 


li- is possible fo categorize the modeling of extended targets by facets into 
those different types: 

i) Modeling by a collection of plane facets of infinite size each of 
which generates a specular type pattern, 

ii) Modeling by a collection of small facets with a uniform pattern 
and size distribution. 

iii) Modeling by facets of finite dimensions with some assumed 
patterns and size distributions. 

2.1 Specular Point Model 

The basic assumption of this model is that the radar return comes only from 
specular points; in other words, it is assumed that we can ignore all the facets which 
are not normal to the direction of propagation. 

This geometrical optics solution has been shown by Kodis [9] to be equivalent 
to the physical optics solution when the principle of stationary phase is used to 
evaluate the Helmholtz integral. In order to show this he performs the integration 
at a typical stationary phase point and averages the result following integration. 
Then the problem resolves itself into one of determining the average number of 
stationary phase points per unit area and the average contribution from each point. 

Kodis applied his method to a rough surface. Then by assuming that the 
return from each specular point has a random phase, he showed that each illuminated 
specular point scatters like a tangent sphere whose radius is the mean of the two 
principal radii of the surface at the specular point. Then, because the radar cross 
section of a perfect conducting sphere of radius r is equal to irr^ (high frequency), 
the return from the target is: 

r = 77 (2-1) 

where r^ and r 2 are the principal radii of curvature of the surface and n^^ is the 
number of illuminated specular points. 
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Barrlck [2] examined fhe geomehrical considerafions ^elevan^ to ihe occurrence 
of stofionary phase points. His analysis is similar to the zero-crossing problem in 
communication theory and his approach follows the same approach that can be found 
in Rice [16], For a two dimensional Gaussian distributed surface Barrick [2] found 
the average points per unit surface area to be 
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(2-2) 


where 


2 

h = mean square height 

^ = correlation distance of the surface 

0 = local incident angle 

ny^= number of specular points per unit area 

From the geometry of the surface, the product of the principal radii of 
curvature can be expressed in terms of the derivatives of the surface of the target, 
C (^ry)r os follows [24]: 

p-3) 

( I 


From the use of the stationary phase technique 


F and 
^x 


are known 


quantities. Hence, the variables in Eq. (2-3) only appear In the denominator, 

Barrick argued that because in getting Eq. (2-1) we use the stationary phase 
technique, it follows that the denominator of Eq. (2-3) is not very small . Hence, 
he assumed that 




( /t- 3 k. -r 5^ ) 


(2-4) 
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Then by averaging ihe denominafor he ob^ained 

^ r, /■*[> = ^ ^ (2-5) 


where 

5, = correlaHon distance 

5 = rms slope 

6 = incident- angle 


Then he subsi-ii-u{-ed<j r^ •'2l> Eq. (2-1) and obtained an expression for 

scafl-ering coefficlenf, cr° given by; 



(2-6) 


Barrick's mehhod of calculafing n^ differs from Longuet— Higgin's [12] mei-hod 
in fhaf the laht-er examines fhe average tata! number of specular points rather than 
the average number per unit area. 

Longuet-Higgin [11] in a special case, assumes that the surface is isotropically 
Gaussian. Then he calculated the joint probability distribution function of the first 
and second derivatives of the surface. By using this probability in a method similar 
to zero crossing technique he calculates the probability of having a specular point 
in an area equal to dx dy. Then by. Integrating this probability he calculates the 
total number of specular points on the surface. His result is as follows; 

■ z 

( -r } (2-7) 


where 


Pt ' 

E( U/ V ) = 

h = 





0 

power spectrum of the surface 

the distance between source and surface 


(2-8) 
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Figure 1 . Average number of specular points 
from Eq. (2~ 2). 



Figure 2, Average ba cks cotter ing cross“section 
per unit area cpo From Eq . (2~ 6). 



Barrick [2] in his invesf■iga^lon assumes fhat fhe surface heighf and ifs second 
derivatives, ( are uncorrelated. Seltzer [20] improves Barrick's derivation 

by first using a correlated distribution function to calculate the probability of having 
a specular point for a given height. Then by averaging this conditional probability 
for a given height distribution he calculates the average number of spec'uldr points per 
unit area. For a bne dimensional surface, he obtained 



where 

= ±h!- 

r 

N = number of specular points 

^ = correlation distance 

2 

h = mean square height 

2.2 Small Scatterers Model 

2,2.1 Spetner's Random~Scatterer Model [22] 

In this model it is assumed that the radar return is composed of the summation 
of the returns from a collection of randomly located scatterers. Then 

r, . f; (2-'o) 


where 


= the density of the scatterer 
= the average RCS of each scatterer 


It is further assumed that the density of the scatterers remains a constant up to 

— 2 

a cutoff wavelength X and then decreases as X for X > X ; 

o o 



x4 CWfo)''^ 


(2-11) 


X > Xc. 
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where 

P = fhe acf-ual density of scotterers 
o ' 

= a constant which is a function of surface properties 
If the effective area of the scatterer is large in comparison to wavelength then: 




( 2 - 12 ) 


where 

A = wavelength 

A = effective area of the scatterer 
s 


But if is small then the RCS will obey the Rayleigh law and will be 





(2-13) 


where C2 is a function of the scatterer size. So the average RCS of a scatterer will 
be 


A" 


(2-14) 


where X is the transition wavelength . 

By substituting Eq, (2-11) through Eq. (2-14) into Eq. 2-10, the results shown 
in Table 1 are obtained. Graphically cr° vs, wavelength is shown in Figure 3. 


2,2.2 Waite's Model [26] 

In this model it is assumed that 1) the target may be represented by a collection 
of scatterers which are located randomly on the target, 2) amplitude of the return from 
each scatterer is a sample of a Gaussian distribution function with zero mean, and 
3) it is further assumed that RCS of scatterers are independent of the frequency of the 
incident wave. 
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TABLE 1 



X<Xq . 

X>Xo 

i 

o 4^p^pA^s 

2 

.o- 2’’°'' 

" ■ 



^ 4jrC.pA^s 



° ■ X^ 



WAVELENGTH 


Figure 3. Wavelength dependence of normalized radar cross~section 
(random scatterer model) from Spetner and Katz, 1960. 
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Based on fhese assumpfions Waire has shown fha^ the distribution of the 
returned power will have the form: 

( ^/,-i 1 



where ^ t 

^ = X 

N = number of independent samples which is equal to one for 
monocromatic case 

th 

a^ = amplitude of the return from the n '' scatterer 
p(w) = probability distribution of the return 


He further showed that: 

N = 


transmitted signal bandwidth 
system resolution bandwidth 


2.3 Facet Model 


2.3.1 Spetner's Model 

In this model it is assumed that the surface is reflective, continuous, and 
with continuous derivatives. It is further assumed that : 

i) The effective mean square slope of the surface is constant if we increase 
the wavelength from zero to some number \ 2 r we increase it 

further, the effective mean square slope starts to get smaller (Figure 4). 
It is further assumed that the density of the facets at normal incidence, 
dsi(O) / varies with the wavelength as follows 


(») =• \ 






>y- ^2. 

(2-16) 


where Cg is a constant close to one. 
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(ii) The radar cross secMon of a large face^ for.verHcai incidence it- 
equal to 



‘rv fli 





V7A2- 


(2-17) 


where is ^he area of fhe facet and is a constant. 
Then, 


O 

- 


4i. Z 


IL 



X 



X 7.’ 



(2-18) 


So from the above assumptions, and when the angular effect is included, 
the return from a target which consists of small and large facets, assuming a Gaussian 
slope distribution for the large scatterers, will be in the form.of 





'Xy'y-T. 


(2-19) 


where 


Cs,C^, 

A 


e 

A 


are only function of surface properties 
average area of large facets 

transition wavelength which is a function of surface 

properties 

look angle 

wavelength 


10 







S = ac^ual mean square slope 

S ^ \2/\ = effecfive mean square slope 

A = con ^ribul■^on from small face is 

The wavelengfh dependence of cP a^ normal incidence is shown in Figure 5. 

2,3.2 Kaizin's Model [7,8] 

This model differs from Spefner's in fhai- he infroduces a paff-ern of finite 
width for the main beam of large facets, thus large facets which are not located 
along the specular direction also contribute to the backscattered return. 

Katzin's first model [7] is primarily concerned with backscattering at large 
incident angles. In fact it is only valid for angles larger than about 85*^. In his 
second mode! the range of validity is broader and it ranges from 0^ to 40° incident 
angles. 

2.3.2. 1 Katzin's First Model: Grazing Angles [7] 

In the first model the facets are divided into two types: 1) large facets, which 
are large compared to the wavelength so that their return is highly directive and can 
be predicted by the physical optics method, and 2) small facets, which are so small 
that the return can be assumed to be proportional to f^ (f is the Incident frequency) 
and the pattern of the return is uniform. In addition, the following assumptions are 
made: 

1) The surface of the target may be considered to be composed of 

a collection of facets of different sizes distributed around a mean 
surface in Z direction. 

2) The facets are randomly located so that the backscatter powers from 
the two types of facets are additive. 

3) Shadowing and multiple reflection effects are neglected. 

4) For a large facet, in small grazing angles, the relation between the 
average radar cross section per unit area and the actual area may be 
derived through the following three assumptions [6] 

(i) ^ plane of facet lies vertically 

(ti) ^ )° plane tilted at an angle 6 

from the vertical position 
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plane Hlfed ah an angle 0 
from fhe verhical and rofahed 
by an angle jd in azimufh 


(in) 


% t.%^ r' 


where o ~ RCS, A = area, and X = wavelengfh. From (ii) and (iii) if follows fhaf 
for facefs for which 0 = cons fan f, and tt > > 0 an average RCS over may be 

obfained, i.e. 

^//l- < (2-20) 


where'o is fhe average RCS for a fixed 0. 

The mafhemafical proof for fhe lasf resulf is given in Aopendix A. 

5) For small facefs if can be assumed fhaf fhe RCS is proporfional fo 
X ^ (Rayleigh region). Then from [19] 



ft 







T 



/v/ 




( 2 - 21 ) 


6 ) 


where 0^ is fhe grazing angle ( 0^ = 'fr/2 - 0 ). 

Because ^/A is an increasing funcfion of A/ X^ for small facefs 
and if Is a decreasing funcfion for fhe large ones, (assumpfions 4 and 
5), Kafzin assumed fhaf we can exfrapolafe cr/A vs. A/ X for 
large and small facefs info fhe infermediafe region wifh fhe infersecfion 
poinf serving as fhe fransifion poinf (Figure 6). 
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_A 


small and large Facets. 



By using the physical optics method, it can be shown, (Appendix A) that 


\r' . ^ 


- 1 / 


/eo (2-22) 

Consider first the horizontally polarized case. As stated earlier for small facets 




1 


I- 


2 

2, 




(2-23) 


The transition point is defined by 




small facets 
from Eq. (2-21) 


large facets 
from Eq. (2-20) 


(2-24) 


where A, is the transition area for horizontal polarization, given by 


(/I.) 


U 



(2-25) 


From the above equation and from the relation between the diameter of a 


circular facet and its area 


D- 



then 


77 D , =• 




r}’^\ 


-z, 4 


X 





— L— i t<ln ^3 

^/r% 4 q 






‘/c 




(2-26) 
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where 



where Is hhe dlamefer for a facef with fransIHon area and g is a number which 
varies from .4 to 1.2 for grazing angles changing from 10° to 30°, 

For vertical polarization we can argue in the same way and get a 
similar equation for 

Figure 6 shows that the normalized radar cross section for facets with area 
equal to A-j is maximum. It follows that the facets for which the diameter is 
comparable to the wavelength have the laigest RCS , So for a target, the bigger the 
number of facets with a diameter comparable to the wavelength the larger the return 
will be (Figure 6). Hence there is a direct relation between the distribution of the 
facet size of a target and the wavelength dependence of the return. 

If the distribution of the slope and the distribution of the area of the target face 
are known we can calculate the backscattering coefficient as follows 



fti 




(2-27) 


where 

and Pg = are respectively, the distribution function of slopes for 
small and large facets 

Ng and = are the distribution functions of facet sizes 

^0'^1'^2~ °^° re;spectively, the areas of the smallest, transition, 

and largest facets 



slope 
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return coniribul'ed by small facets 
return contributed by large facets 
average RCS of a small facet from Eq. (2“21), 
average RCS of a large facet from Eq. (2-20) 

2. 3. 2. 2 Katz in's Second Model: High Depression Angle [7] 

For near vertical incidence the radar cross section of large facets computed 
by using the physical optics method is 




( 


(< vw'F 



(2-28) 


where 


y = si n 6 - Z cos 6 

y 

5 = -T cos e 
cos X = “COS 0 - Z sin 0 

y 


and 

jt, w = are dimensions of the facet 
A = the area and is equal to Jl x w 
Z^,Zy = are, respectively, slopes of the facet in x and y directions 
\ = wavelength 

For this model the RCS of the small facets, for small incidence angles, 
is assumed to be 


{.^5 ) 





(2-29) 


With the RCS for the two types of facets known, the total RCS can be obtained 
by averaging each individual RCS over its size and slope distributions and then sum 
them together, i .e . 

rs . (2-30) 
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where 


fi" A) 10 ) A 

K,vjV ' 


-<* A, 


^ Y.Jlh ^ ^ 


Art, A,, Art = are respecfively, hhe smallesf, fransiHon, and fhe largest 
are in facet distribution 

P(Z^, Zy) = distribution function of slope 

N(A) = distribution function of area 

In order to calculate the above equation Katzin assumes that 

(i) kiA/ ^ ^ 


(ii) <;o^ >1 

2 2 I 

(iii) 2K-^<3 o<s»' I ^ ^ 
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where 

K= ^-^ 1 % 

w = width of facet, x direction 
A = length of facet, y direction 
6 = incident angle 

Cy, = are respectively, the standard deviation of the slope along x 
and y directions 
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3.0 PROP.QSED FACET MODEL 


Data taken from disks of different sizes shows that the variation of the RCS 
as a function of the disks dimensions has the shape shown by the solid line in Figure 
7 below; where Ka = 2Tr/x- x radius of the disk [19]. It is apparent from the solid 
curve in this figure that there are at least three recognizable regions, and hence 
it appears appropriate to separate the size of facets into the three categories rather 
than two. Approximately, the following three types of facets may be defined; 

(i) Large facets for which Ka>ir and the return can be determined by the 
physical optics method. 2^^ 

(ii) Medium size facets for which tt> Ka> 5 • For ^^bis region the return 

is almost constant and does not change more than 3 dB for the whole 

range. It may be assumed that the value of the RCS for the whole range 
. ... *5/4 

vanes with A . _ 

2tt 

(iii) Small size facets for which 5 > Ka. For this case the small perturbation 
method can be used to evaluate the RCS. 

By dividing the facets, according to their sizes, to three instead of two 
categories, we overcome the problem of extrapolating the curves for large and small 
facets to intermediate sizes. This should give us an improved result over the two 
scale model; since the assumptions on facet sizes are less severe. To see this note 
that in reference 19; 
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TT 


K a. 




(3-1) 


Hence, for 




(3-2) 




o 



(3-3) 


It follows that Ka= 1. The dashed line in Figure 7 shows that Ka= 1 belongs to the 
small region and by this extrapolation not only the entire intermediate region but 
part of the small region are grouped into the large size region. 
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Figure 7. RCS of a disk as a function of deminsion/wovelength . 
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Katzin assumes fhaf the probability distribution function of the area is in the 


form of 






(3-4) 


where and m are constants and P(A) is the area probability distribution function. 

In order to overcome the problem of determining P(A), we start from the 
dimensions of the rectangular facets instead of their areas. Then the three regions 
will be as follows: 

(i) Small facets for which dimensions are smaller than 2X/5. 

(ii) Medium facets for which dimensions are bigger than 2\/5 but smaller 
than X . 

(iii) Large facets for which dimensions are bigger than X . 

In order to compute the return we have to know the distribution of the dimensions 
of the facets. This distribution can easily be calculated by using the distribution 
function for height and slope of the target; Thus, an independent assumption about the 
facet area distribution can be avoided. 

For the facets with dimensions less the 2^5 we' suggest to use the small per- 
turbation method. In this case it can be shown that (Appendix B) 

- (3-5) 

where 

2 

h mean square height 
Z length of facet (max.) 

We should add that the answer from this model has the local incident angle and if we 
want to consider it as a superposition of the small facets over the large ones, we 
should average it over the variations of the angle due to the variation of the slopes 
of the large facets. 

In what follows we calculate the RCS of the large, medium, and small facets. 

It is assumed that the distribution of the facet size is in the form of 
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If is also assumed thaf fhe slope has a Gaussian disfribufion with zero mean. 

3,1 Large Facets 

From physical optics the radar cross-section of a plate with a length of "a" 
and a width of "b" is equal to 

Now if we assume that the direction of propagation is in YZ plane and has 
an angle 0 with the 2 direction, and the plate has a slope Z in x direction and Z 
in the y direction than from [ 7 ] the cross section of an arbitrarily oriented rectangular 
plate is 
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rt 


(3-7) 

Kb \ 

Kb / 


where 

k' -_?2L 
^ X 

a,b, are respectively length and width of the plate 


z ,z 

x' y 


are respectively the slopes in x and y directions 


In order to get contribution of the large facets to the return power we should 

average (0) over all possible slopes and all possible size. We assume that the 

slope distribution in both x and y directions are Gaussian with zero mean and 

standard deviation of s and s . Then 

X y 
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where 


ffL(«) 


contribution of the large facets 


P(a)/ P(b) 


are respectively the probability distribution 
functions of the length and width of the plate 


'1' 


are respectively the smallest length and width for a 
plate to satisfy the physical optics requirement 


As stated earlier we assume that the distribution of the size of the facet is a 
Gaussian function, i.e. 
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By assuming that 

i) 

H) A. = 

The contribution of large facets will be 


(3-10) 
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3-11) 


where 


B = the size of the largest facet on the target 

The integral over and have the general form of 

^ / JUm. PCc-1- ^ (3-12) 

V p(Ci-6^] ^ 

-tA 
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where the parameters are as follows; 


i) integral over Z 


X 


ii) integral over Z 


5 - ^ 0 




(3-13) 


/i = 




p» K- h 


(3~14) 


Appendix D shows that the answer to the integral in Eq. (3-12) is 
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(3-15) 


where 



(3-16) 


Z = X + iy 
Z* = X - iy 
X =PB/A 
y =AC/& 
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From Appendix D 
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(3-17) 


Now if we consider 1-he lost integral of Eq. (3~1 1) we have 
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We see l-hat the first integral in the right hand side has the form of Eq. (3-12), the 

second is zero, and the third one has the form of Eq. (3-17). So in order to calculate 

this integral we first calculate it from Eq. (3-15) then differentiate it vs. A and 

get _ . 

9A2 

For the case when the parameters x of the error function is large i.e. when 
we are dealing with very large facets in comparison with wavelength, and these are 
the facets which contribute the most, we can assume that the error function is equal 
to unity and this will make the computation much easier. For this case 


(i) Integral over 
From Eq. (3“1 1) 


integral over ^ 

where I is defined in Eq. (3-12). Then by using the result of Eq. (3-15) and by 
assuming that the parameters are as indicated in Eq. (3-13) 



(3-19) 


where 

X = Ka V2~ S COS0 
X 

K = wavenumber 
a = length of facet 

S= standard deviation of slope in x direction 
6 = incidence angle 

2 

“X 

then from here because it is assumed that x>2 so erf(2) ^ 1 ; e 0, 













(3-20) 
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Now if we further assume that 

/2“ S K cos6»l, i.e. facet size is large and the standard deviation of (3-21) 
X a 

the slope is not very small. 


integral over Z^= 


(3-22) 


ii) Integral over 

From Eq. (3-18) and by using the results of Appendix D as indicated in Eq. (3-17) 






(3-23) 


where I is defined in Eq. (3-12) and its parameters are given in Eq. (3-14). Then 
by using the results of Eq. (3-15) with the same procedure as indicated in getting 
Eq. (3-22) 
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Now if we assume that 

1) x>>! 

2) x> y 


(3-^25) 

(3-26) 


i.e. consider large facets only near specular direction, then 
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By substituting these results in Eq. (3-11) 


6 



. (3-30) 
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assuming that B ->- «> , i.e. the size of the largest facet is very large. 
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(3-31) 


This result, as it is expected, is in general form of the result that one would get 
by using the geometric optics method. The only difference is that we have (!-► ® j 

rather than ( /-+ . The reason for this is the assumption that is made in 

Appendix D, i ,e . 
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^ V p(c+8v>) 


(3-32) 


which is not true for the extreme case when B and the angle of incidence is 
not equal to zero. 

It should be noted that in comparing the theory with the measured data we 
do not use the assumptions indicated in Eqs, (3-25), (3-26) and (3-31). For this 
purpose we use the result of Eqs. (3- 15) and calculate it numerically or digital 
computer. 

3.2 Medium Size Facets 

The facets whose dimensions are less than X but larger than 2X/5 are con- 
sidered to be medium size facets. 

There is no exact easy solution to the RCS of the medium size facets. Figures 

8 and 9 show that for small and large facets, qs it is expected, the normalized 

o 

RCS versus area varies, respectively, with A and A (A Is the area). We also see 
that the slope of the curve for medium size is smaller than the slope for both large 
and small size, and it varies almost as a . Figure gives us the variation 

of the return versus the local angle of incidence. This variation is in the form of 
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Figure 8. Measured cr° for small and Figure 9, Measured < 7 ° for small, medium, 
medium facets. and large facets. 



0 15 3.0 45 60 75 90 0 15 30 45 60 75 90 

Incidence Angle in Degrees Incidence Angle in Degrees 


Figure 10. Measured W polarization backscatter cross“section versus the angle of 

incidence for perfectly conducting rectangular planes with heights of 2 %q 
and lengths of 0.5 and 
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So in view of fhe meosured dafa shown in Figures 8,9, and 10 we assume 
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where 

K =2ttA 

X = wavelength 
a = width of the facet 
b = length of the facet 
9' = local incidence angle 


If We write Eq . (3-33) in terms of incidence angle 
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where S is the slope 


(3-34) 


Because we assume that the slope of the medium size facets are not large and because 
the slope of the return versus angle is small for the medium size facets, then s 0. 
Hence Eq. (4-34) will be: 
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(3-35) 
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In order to calculate the contribution of the medium size facets w'e have 
to average Eq. (3- 35) over a and b. Then 


> 



5 


where 

K =27tA 

A = wavelength 
a = width of the facet 
b = length of the facet 

p(a) = distribution function of the width of the facet which is assumed 
to be Gaussian 

p(b) = distribution function of the length of the facet which is assumed 
to be Gaussian 

0^(0) = contribution of the medium facets to RCS 

In order to compare the result of the theory with the data we have to take 
into consideration the reflection coefficient of the facets. The assymptions^for 
calculating the reflection coefficient has been explained in Section 5.1 . 


3.3 Small Facets 

In order to calculate the contribution of small facets i.e. facets whose 
length is less than 2A/5 we use the perturbation method as first formulated by Rice 
[17] and then further developed by Peak [14] . 

In this technique unlike the tangent-plane approximation there is no require- 
ment indicating that the radius of curvature at any point must be large In comparison 
with the wavelength. So we can easily use it for the return from small facets of the 
natural surfaces. 

The requirement of this model are as follows [19]: 


i) h/^ < ( 


^ i.e. roughness height is small 
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i.e. surface slopes are relafively small 



i.e, roughness is Isotropic 


In. the previous section we showed that the first requirement is correct for 
small facets and by assuming that the second and third requirements are also correct, 
the radar cross section will be [19] 


^ B- HI X (3” 37) 

where 

X = wavenumber 
h = r.m.s. height 
6 = angle of incidence 


In this equation a, for horizontal, cross, and vertical polarization are given, 
respectively, by 
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(3-38) 
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and I is defined as. 



j V [2 K- Xc(?' r ) 
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(3-40) 


where J|^( ) is the zero-order Bessel function and /^(r) is the surface-height correlation 
coefficient for an isotropically rough surface and is defined as 


/’CO. 




(3-41) 


where 

yr, J iy % ) 


(3-42) 


For the case that the surface is not isotropic the result for Eq. (3-40) is given in 
[ 5 ] , Now if we assume that the correlation coefficient has a Gaussian form, i.e. 





(3-43) 


then from Eq. (3-40) [19] 
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(3-44) 


where Z is the correlation distance. 
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By substituting Eq. (3“44) in Eq, (3-37) and assuming - 1 
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4.0 PROPERTIES OF THE SYSTEM 


In fhis chapfer we explain Ihe principle and ^he consIrucHon of fhe large 
fime-bandwidi-h FM“CW radar which is used fo colled dala from ihe bare ground. 

Sedion 4.1.1 is devofed fo briefly state the principle of the radar. In 
Section 4.1 .2 we have calculated the frequency spectrum of the radar and have 
shown that it is uniformly spread over the sweeping bandwidth. This result has been 
used in Section 4.2 to calculate the variance of the return power. 

Section 4.2 examines the effect of the frequency averaging. We first show 
the predicted reduction of variance due to frequency averaging by assuming a uniform 
scatterer model [26], then we extend this result to take into consideration the pene- 
tration effect. Then we calculate the reduction in the variance for the case that the 
target is assumed to consist of a collection of facets. 

Section 4.2.3 compares the theoretical results obtained in Sections 4.2.1 
and 4.2.2 to the measured data. It is shown that by considering penetration we can 
explain the reduction in the variance of the backscatter data from alfalfa. We also 
show that the facet model can explain the difference between the measured curve from 
the bare ground and the curve predicted by assuming that the ground consists of small 
scatterers . 



4.1.1 Radar 

The radar used in this investigation is a FM-CW radar which has a very large 
time-bandwidth product. A FM-CW radar transmits a modulated signal as shown in 
Figure 11 . Then from [21] the IF frequency after the mixer is equal to 



C 




(4-1) 


where 

Fjp = intermediate frequency 
R = • range 
c = speed of light 
F = RF frequency deviation 
F|^ = frequency of modulation 
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A block diagram for ihe radar used Is given in Figure 12 and ifs characferisMcs 
are as follows [3]: 


TYPE 

FM-CW 

Modulating Waveform 

Triangular 

Frequency 

8-18 GHz 

FM sweep: Af 

400 MHz 

Transmitter Power 

10 dBm (10 mW) 

Intermediate Frequency 

60 kHz 

IF Bandwidth 

3.58 kHz 

Antennas 


Height above ground 

26 m 

Reflector Diameter 

61 cm 

Feeds 

Cavity backed. 


log“periodic 


Frequency 

(GHz) 

Calculated 

Antenna 

Gain 

(dB) 

Effective Beamwidths of Product Patterns 
(Degrees) 

Az El 

8 

31.2 

2794“ 

3.43 

10 

33.0 

3.07 ■ 

3.24 

12 

34.6 

2.42 

2.38 

14 

35.9 

2.35 

2.34 

16 

37.1 

1.65 

1.46 

18 . 

38.1 

2.02 

3.20 


4.1.2 Frequency Speclrum of ihe Signal 

In ^his secHon we would like fo show ^ha^ ihe frequency speclrum of fhe refurn 
signal of an FM“CW radar wilh which the data has been collected is flat over the 
sweep band. If the frequency of the wave varies with time as shown in Figure 13, 
then the wave form will be: 


RECEIVE 

ANTENNA 



• Figure 12. Basic block diagram of S’- 18 GHz radar specfromefer. 







Figure 13. Frequency of the transmitted signal as a function of time. 
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Because s(f') is a periodic function we can wri^e 
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First- we should calculate F^. By substituting Eq. (4“2) in Eq, (4-5) we get 
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By using the relation 
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Eq . (4-8) becomes: 
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From ^he Fresnel integral properMes 
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and because fhe system used has a very large Mme-bandwidth produci", i.e. 
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It can easily be shown [4] that from Eq. (4-10) the signal has spectrum only around 
the carrier frequency, and its magnitude is 
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So from Eq . (4-7) 
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M =-«0 


This spectrum is shown in Figure 14. 
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Figure 14. Frequency spectrum of the radar. 
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4.2 Frequency Averaging 

It is a well known fact that radar return from irregular targets may fluctuate 
due to insufficient averaging. In order to overcome this problem, over the years 
many techniques have been used, and frequency averaging is one of them. Waite 
[26] has shown that if we average the return of the target for different frequencies 
the standard deviation of the average will be much smaller than the standard 
deviation of the data measured in each frequency. For doing this, Waite assumed 
that the target is composed of small scatterers for which the pattern is assumed to 
be uniform. 

In Section 4.2.1 we shall first explain the work done by Waite then we 
shall expand his work by taking into consideration the penetration of the wave into 
the target. 

In Section 4.2.2 we assume that the target is composed of the facets rather 
than small uniform scatterers. We have shown that if the return from different 
facets ore comparable than the reduction of standard deviation of the return due to 
frequency averaging will be the same as predicted by Waite. But if the return from 
different facets are not comparable, then we have to have some information, about 
the facets whose contribution to the return is very large, in order to be able to 
determine the reduction of the standard deviation of the return . 

In Section 4.2.3 we compare the theoretical results obtained in this chapter 
by the actual measurement and show that alfalfa is a good example of a target composed 
of layers of small scatterers. We then use the results of the facet model to explain the 
number of independent samples in the RCS measurement of the bare ground. 

4.2.1 Uniform Scatterer Model 

4 2. 1.1 Uniform Scatterer Model with no Penetration [26] 

In this model it is assumed that the target may be represented by a random 
collection of discrete independent scatterers. Each scatterer has a uniform pattern 
and its reradiation amplitude is a sample of a Gaussian random process with zero 
mean and standard deviation. 
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As Figure 15 shows fhe illuminal'ed area of the ^a^ge^ is a funci*ion of D and 3 , 
fhen ihe refurned signal will be in ^he form of 




where 


s(f) 

r 

n 

a 

n 
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roundfrip lime delay of n scafferer and is equal to 
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The total return from the illuminated area is: 
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From the characteristic of delta function we can write the transmitted signal 
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Hence, the return signal may be written as 
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where the symbol ® denotes convolution of two function [10 ] . 


(4-16) 


(4-17) 


(4-18) 
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(4-20) 
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The fargef effeci’ is only inlroduced in 1‘he second farm of ^his convolu^ion so if we 
define c(l-) as 


C[}) ^ ^^(.+■'■^'1) 


(4-21) 


then c(t) will be the target impulse response. 

In order to obtain the reduction in the standard deviation of the return signal 
due to frequency averaging, consider the Fourier transform of c(t). 
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From Eq. (4-22) the autocorrelation of the response of the target in the frequency 
domain can be calculated as: 
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From fhe requiremeni- of the model that the scatterers are independent samples of a 
Gaussian process with zero mean we obtain: 


Hence 
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By the assumption that a^ and t^ 


are independent 
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By assuming a uniform distribution for time delay i.e. 


^ 

1 _ 
T 

r* --I 


1 e 





^ ^ - 1 ^ 2 . 






X 


y. 


s- 

T 


'T. 


— * 


(4-28) 


(4-29) 


51 




*[ 

which indica^es fhat for fwo frequencies -:j|s- Hz oporf fhe fargef effecf will be 
uncorrelafed. 

In whaf follows we will obfain fhe mean and fhe variance of ihe refurned 
power and by using fhe correlation of the target in frequency domain from Eq. (4“30), 
we will calculate the reduction in the variance of returned power due to the frequency 
averaging . 

Variance of the Returned Power 

As' stated earlier in Eqs. (4-20) and (4-21) 

C+) 51+) ® CC+J (4-32) 


Hence, the Fourier transform of S^(t) is 


(f)= S (f } cif) 

then the received power will be equal to 
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Now by using this relation we will calculate the variance of the return power, 
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Now we calculate each term separately 
i) Mean 
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ii) Expected value of w 
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In order fo calcula|-e ihis .Integral we use the procedure as. stated in reference [10] 

t .e . 




( 4 - 38 ) 
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Now by using the results of Eqs. (4-36) and (4-40) we calculate the variance of the 
returned power, 
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Because, as shown in Eqs. (4“46) and (4-47) 
is no^ a funcMon of frequency 
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In whaf follows we will calculafe fhis infegra! by using fhe assumpfion ihal Ihe t'orget 
is composed of small scaUerers with uniform paffern. In reference [26] if is sfafed 
fhaf in this case 
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This result’ is correct' If fhe following assumptions are satisfied; 


i) Number of the scatterers contributing to the return is large, 
ii} fT, where f is the frequency and T is the time span of the target 
is very large. 

Then by using Eqs. (4-23) and (4-30) we obtain 
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and 


c(f) c" ^ ' 


Because both of these results are independent of the frequency we can use 
Eq , (4“44) to obtain the variance of the returned power 
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where p(u) is defined in Eq. (4-43). 
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Now in order ho calculafe hhe improvemenh due ho hhe frequency averaging 
we compare hhe shandard deviahion of hhe rehurn when hhe frequency is sweph over 
fm wihh hhe case when hhe sweep is zero. 
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In Sechion 4.1.2 ih was shown hhah hhe signal has a frequency spechrum as 
shown in Figure 16^ where hhe power of hhe incidenh signal is assumed ho be one 
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fhen from Eq. (4-43) 
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Figure 16. Power spectrum of transmitted wave. 
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or 


(4-52) 







In order fo calcula^e Eq. (4“49), we firsf calculofe fhe standard deviation 
of the return signal by substituting Eq . (4-52) in Eq, (4-48) 





(4-53) 


For the case, when there is no frequency averaging fm 0, i.e. 
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Therefore, fromEcj. (4-49) 
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So 



I y^^M* T" 



(4-56) 


Where I is the reduction ratio of the variance that we get when we perform frequency 
averaging (Figure 17 ). Note that from Eq. (4-54) and Eq. (4-44) 




no fr/e. 




(4-57) 
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Variance Reduction Ratio, I, Due to Frequency Averaging 



Figure 17. . Variance reduction with rectangular spectrum 
panchromatic illumination . From Eq.4~ 56. 
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so Eq. (4”55) can be w^i^^en as 
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1 /9/e, 




(4-58) 


So by rhe assumpfion of small sca^^ere^s we can calculate fhe average y and hhe 
deviaHon a of hhe dafa teken by a FM-CW radar, and fhen fhe variance reducHon 
due fo ihe frequency averaging will be 


r- 



(4-59) 


4. 2. 1.2 Uniform Scaijerer Model wifh PenefraHon 

In ibis secHon if is assumed fhaf fhe fargef can be represented by independent- - 
scatterers located in layers. The pattern for each scatter is uniform and its return 
is a sample of a Gaussian random process wifh zero mean and crstandard deviation. 

The procedure for calculating the variance reduction due to frequency averag- 
ing is very close to the method used in Section 4, 2. 1.1. We first calculate the 
correlation of the target return in frequency domain then use this correlation to 
calculate the variance of the return signal. 

From Figure 18a 
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where 


e 

a 

h. 

4 , 



incidence angie 

beamwidfh of anienna 

heighf of fhe source from fhe i*^^ layer 

smallest- distance between the source and the illuminated area 
of the i*^^ layer 

largest distance between the source and the illuminated area 
of the i^^ layer 


then 

T,,= -ih; 

2.kt 

C 


where 

C = speed of a> light 

T.-. = smallest time required for travel, a roundtrip distance, 

il 

between the source and the i layer 

T,„ = largest time required for travel, a roundtrip distance, 

tfl 

between the source and the i layer 
As stated in Eq. (4-1) the return from the n*^^ scatterer of the i^^ layer is 


(4-62) 

(4-63) 
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Then the total return will be 




(4-65) 


64 



w?^h bhe same procedure fhah is used, to get Eq. (4“22) 
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then by the definition of Eq. (4-23) 
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then by the assumption 
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we have 
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RC>^] - E I 2T ^ 2” e -r • 



(4-71) 


where 


R( v) is fhe correiafion funcfion of fhe refurn in frequency domain and 
/ e is ^he contribufion of i layer 


By assuming fhaf ihe Kme delay of fhe rel-urn from each sca^^erer is a sample of a 
uniform disfribuHon, i.e. 


f u- ) = 



(4-72) 


ofherwise (4~73) 


and by assuming ^ha}■ fhe return from the scatterers of each layer are comparable in size 
before attenuation through layers above them is accounted for, we obtain from Eq. 
(4-71) 
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by the same procedure as stated in Eq. (4“29) 
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so in view of Eqs. (4-60) and (4-61) we can write 
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where 
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ih 

The incident power that reaches the i layer has already been attenuated by the 
first (i“l) layers. The power which is scattered from a scatterer in i ^ layer is 
also attenuated by the first (i”l) layers before it reaches the receiving antenna. 

So assuming that the average RCS of all layers are equal the contribution of the 
lower layers will be much smaller. From Figure 18b it can be easily seen that the 
distance in which the power gets attenuated is ( hi - h , ) q' where 6' is the 
refraction angle and from the SnelPs law is equal to 


/ 


e = 





(4-78) 


where 

6 = angle of incidence 

= relative dielectric constant of the target 

So if we assume that the attenuation function is exponential, then 



Chi-hi) 
S' c<r^9' 


(4-79) 


where 6 is the skin depth and a function of the dielectric properties of the target. 
By substituting Eq, (4-79) in Eq. (4-76), we obtain 
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By subs^it•uf^ng Eq. (4-77) Info Eq. (4f-80) 



Because the distance between source and surface is large, the skin depth is very 
small, and the refraction angle as defined in Eq. (4-78) is small, so the area of 
illumination for different layers is almost constant. Taking h. to be equal to h-j 
except in the exponent, we obtain 






The above sum can be approximated by the following integral 
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Where 


N — number of fhe layers 


(4-85) 


N of J'he variance of all scatferers locafed (4-86) 

^ in all layers . , . 

In order fo calculafe fhe standard deviation of the average return signal we substitute 
Eq. (4-84) and Eq. (4-52) into Eq, (4~48). We obtain 
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(4-87) 


where R( v) is defined in Eq. (4-^4), 

The reduction of the standard deviation due to frequency averaging from 
Eq. (4-49) will then be:. 
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From Eq. (4-84) 
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? ^<5-^ 

So Eq, (4-88) will be 
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where 


T 






is fhe Hme span of fhe fargef , Now if we lef fhe depfh of fhe fargef go fo infinity, 

i.e. 

h — *■ oo 


max 
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Then Eq. (4-90) will be 



where 

T = ihe Hme span of fhe target 

6 = skin depth 

c = speed of light 

9 = angle of incidence 

6 = antenna beamwidfh 

V = sweep frequency 

f = sweep band of the signal 

m 

I = improvement factor 

9 = refraction angle 

Because the requirements leading to Eq. (4-59) are correct for the layer assumption 
so for the measured data the reduction ratio of variance due to frequency averaging 



(4-92) 


where 

CT = variance of the measured data 
y = mean of the measured data 
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4,2.2 Number of Inc^ependen^ Samples Based on Facefr Model 

In this section we assume that the target consists of a collection of facets 
of different sizes, and the return from different facets are independent. 

In the previous section it was shown that the target effect can be represented 
in the form of: 

C(f) = T 


-3 

C 


(4-22) 


The mean arid the variance of the return power can be expressed in terms of c(f) as 
follows: 





(4-54) 







It was also shown for the uniform scatterer model, that in the absence of 
frequency averaging, the variance of the return power is equal to its mean, i.e. 


S [ 1 ''J 




(4-93) 


Based on this assumption for calculating the reduction factor of the variance 
of the return power due to frequency averaging, instead of conducting two experiments, 
one with frequency averaging and the second at a single frequency, and calculating 
the variance of each set of the data, we only need to know the mean and the variance 
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of fhe firsh sef. The reason is thar because we are sweeping over a narrow band, 
compared wifh hhe center frequency, so the mean with frequency averaging will 
be equal to the mean that one would get from the second experiment. On the other 
hand from Eq, (4-93) the mean square is equal to the variance for the single frequency 
case. Hence based on the uniform scatterer model, if we conduct one experiment with 
frequency-overaging-'then the reduction factor'of the variance can be obtained by 
dividing the mean square of the data by its variance. That is: 






(4-94) 


'wir 


where N is the number of independent samples . 

The experimental data shown in Section 4,3 are from reference [3] and are 
calculated based on this assumption . 

In this section we will first show that when the target is assumed consists of 
a collection of facets, the assumption of the equality of the mean and the variance 
of the return power in the single frequency case does not hold— this is especially 
true for the case that the product of the sweep band and the time span of the target, fmT, 
is small, i.e. at small incidence angles. So we will conclude that the calculated 
number of independent samples as shown in Figure 21 are not correct; because they 
are based on Eq . (4-94). ' 

In the second part we will calculate the number of independent samples 
based on the facet model assumption and show that the results for small incidence angles 
are almost the same as the results predicted by the uniform scatterer model assumption. 
For large incidence angles the number of the independent samples predicted by the 
facet model is smaller than the number predicted by the uniform scatterer model and 
is closer to the number calculated based on the experimental data Figure 21 . 

4.2.2. 1 Mean and Variance of the Return 

In what follows we will first calculate E[ I c(f) I I c(f-\; ) | ] then by putting 
V = 0, we will obtain E[ | c(f)| . We will then calculate the mean of the return 

power, i.e. E[ | c(f)| and will show that, unlike the uniform scatterer model, the 
mean square is larger than the variance, i .e. 
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(4-95) 
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from Eq. (4-22) 
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Because if is assumed fhaf and f^ are independent and f^ in uniformly distributed, 
we get 
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Now we will calculat'e E[ |c(f) 1 ] E[ |c(f“v)l ] and by subs^itu^Ing v = 0 we 
get the mean square. 
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Now we want to prove that in general mean square is larger than the variance, i.e. 


(mean)^ - variance > 0 
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If we substitufe \) = 0 in Eqs. (4-97) and (4-98), and assuming xt’haf N is large 
and 
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As will be shown < a > / < «„ > /< °n ^ ^ °n ^ comparable 

magnitude, so for the small incidence angles where the time span of the target, T, 

is very small fT will not be large and consequently j ^ 

is not very small, and because N is large, the first term in Eq. (4-101) will 

dominate . So 
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In other words the mean is larger than the variance. 
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4. 2. 2. 2 Number of Independeni' Samples 

We want to calculate the reduction in the variance from Eq. {4"44). From 
Eqs. (4“97) and (4~98) the integrant in Eq. (4-44) will be: 
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From Eqs. (4-44) and from the definition of the number of the independent samples 
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We obtain 
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(4-104) 

where p(v ) is defined in Eq, (4“52). When the product of sweep band and time 
span of the illuminated area is small i.e. v T« 1, the answer to both integrals 
will be one, so I, as expected, will be equal to unity. But when the product is 
not small and since f >> v , fT will be large and j ^ 0, so I will 

be equal to 
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where I' is the number of independent samples predicted by uniform scdtterer 
assumption (see Section 4.2.1). If we substitute N < > by summation in 

Eq. (4-106) we obtain 
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(4-107) 
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A/ 

It should be mentioned thot the summation 2T for ail facets-” large, medium/ 

n= 1 

and small — so in order to calculate the second term in the denominator we break it and 
group each size separately 
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(4-109) 


where 

n^^ = number of large facets 

n„ = number of medium facets 

m 

n„ = number of small facets 

s 

N “ n„ + n + n 

£ m s 

Now we calculate N, n „ , n , and n . 

X* m s 2^ 

Time delay due to each facet with a length L from Figure 19 is — — sin 0^ 

where 6^ is the local incidence angle. The average value of this delay will be: 
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where 

s = slope 

p(s) = probability disfrribufion of slope 

p(L) = probability distribution of facet length 

Since the -time delay is-equal to-the average tirfie delay times the total number 
of facets, we get 
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where, as before, it is assumed that the targets whose size is, smaller than 2 X/5 
are considered to be small, larger than X are considered to be large, and between 
2X/5 and X are considered to be medium. 

In what follows we calculate E[a^^] and E[a^^] for large, medium, and small 
size facets and substitute the results in Eq. (4-108) to calculate the number of the 
independent samples, 

i) For large facets note that 
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Then assuming a Gaussian distribution for facet length, we obtain 
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For relatively large incidence angles 
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ii) Medium size facets 


’ KL yLi<y 


(4-125) 




(T/v^ 


, yii^ KLJ^ 

i, ^ 


<lL 


S 


S' 8 ^ 


.^ , • -j/ 

^ p iff ,, 

— 

L- 

r 

/ .17^ . z (4-126) 

" ^ ^ r 

^ — xJ^ J - ^(-rKX^^j 

( '7 ^ L 


89 



i 


^ S 


S' JTn 

e. 


where 




I^L ^ j 


e a 


'2^ 


' g. , 

3 / 


?k ^ •"^-+ (4-127) 
<■ X. 1 .«K >.^'<^1 


5t 




(4-128) 


and 




!\J 


- I 


(4-129) 


iii) Small fdcels 
From reference [19] 
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where 
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(4-133) 


(4-134) 


Now if we subs^i^u^e fhe results of Eqs. (4-131), (4-132), (4-127), (4-126), 

. (4-118), (4-117), (4-115), (4-114), and (4-113) into Eqs. (4-108) and (4-109) 
we obtain ^ , and by 

substituting these results in Eq , (4-106) we calculate the number of the independent 
samples . 

In Figure 21 we have shown one typical result of this calculation. 

4.2.3 Comparion of the Results with Measured Data 

In this section we compare the results obtained in Section 4.2 with 
measured' data . 

An experiment was conducted to measure the effect of frequency averaging on 
the reduction of the variance of the returned power [ 3] . The results of this experiment 
for alfalfa and bare ground are shown, in Figures 20 and 21 . 

4.2.3. 1 Alfalfa 

At incidence angles smaller than 60° alfalfa can be considered a layered 
medium consisting of uniform scatterers. In Figure 20 we compare the theoretical results 
obtained assuming no penetration, Eq. (4-56), and the case that we have some 
penetration, Eq. (4-91), with the result of the experiment. It can-be seen that for 
incidence angles smaller than 60° the penetration model can very well predict the 
experimental result. For incidence angles larger than 60° the alfalfa behaves like a 
collection of facets, so the number of the independent samples should be calculated 
based on the facet model assumption, Eq. (4-107). 

4. 2. 3. 2 Bare Ground 

In Figure- 21 we compare the result obtained based on Waite's [26] derivation 
and the result obtained based on the facet model theory, Eq. (4-107). It is clear that 
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Figure 20, Number of independent samples for alfalfa 
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at small incidence angles, the facet model does not agree with the calculated number 
of independent samples based on the experimental data. This may be due to following 
reasons: 

i) At small incidence angles, as was shown in Section 4.2.2, the 
return power does not have a Chi-square distribution.. Hence the 
number of the independent samples calculated based on the 
experimental data. Figure 21, are not correct. 

ii) The theory of the reduction of the variance of the return power based on 
the facet assumption is not complete and needs more refinement. The 
derivations in Section 4.2.2 consider only the two dimensional case, 
although the actual target consists of three dimensional facets. This 
deficiency can introduce errors, particularly at small incidence angles. 
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5.0 RESULTS 


In this chaper we will compare the theoretical results obtained in Chapter 3 
with the experimental data obtained from the bare ground. . 

Because in general the target is not perfectly conducting, in Section 5.1 we 
will calculate the reflection coefficient of the target as a function of its dielectric 
constants. 

In Sections 5,2 and 5.3 we will compare the prediction of the model developed 
in Chapter 3 and Katzin's model [8], with the measured data. It will be shown that 
the new model predicts results better than Katzin's model. 

5.1 Reflection Coefficient 

5.1.1 Large Facets 

From reference [23] the reflected field from a piqne boundary between two 
media is 
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At K, 
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where 


k2 = /<,6, 5 

1^2 ^ e, /-iT, ») 

.^•j = Complex permittivity of the first medium 

^2 = Complex permittivity of the second medium 
= Complex permeability of the first medium 
^2 = Complex permeability of the second medium 
For our case = ^/So from Eq. (5~1) 
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and from Eq. (5“2) 
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In order to calculate the reflection coefficient for the ground we use the 
moisture and the type of the soil to calculate e' and €" from Figures 22 [25] 
and then use these results in Eqs, (5-3) and (5-4). The calculated reflection coefficient 
for HH and W polarization are shown in Figure 23. 
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DIELECTRIC CONSTANT 


REAL PART 

A ABILENE CLAY LOAM, 10.6 GHz (WIEBE, 1971) 

□ AMARILLO FINE SANDY LOAM, 10.6 GHz (WIEBE, 1971) 
70r o RICHFIELD SILT LOAM, 9.4 GHz (LUNDIEN, 1966) 

V LOAM, 9.5 GHz (LESCHANSKII et al., 1971) 

• WATER, 9.3 GHz, 20° C (PARIS, 1969) 



SOIL MOISTURE (PERCENT BY WEIGHT) 


Figure 22. Measured dielectric constant data of loamy soil as a function of moisture 
content by weight around 10 GHz. Solid curves were drawn to fit the 
data points and the broken curves were extrapolated [25] . 
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LOSS TANGENT 




Reflection Coefficient for Large Size Facets for Different Moisture Contents 



Figure 23. Reflection coefficient for large size facets. 
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For the case of bare ground we can always assume that = 1 . Then if we define an 
effective reflection coefficient, we obtain 









In Figure 25 [R HH| and |R VV( are shown as a function of incidence angle 
s s 

and soil moisture content. 


5.1.3, Medium Size Facets 

For the medium size facets we assume that the reflection coefficient is the 
average of the reflection coefficient of the small and the large size facets, i.e. 


Rrt ItH 
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and 

v'v' j 

9 9 

Figure 24 shows |R HH| and jR VV] versus angle of incidence for different 
M M 

soil moistures. 

5,2 Comparison of the Results of the Theory With Data 

In the following pages we will compare some samples of the measured data 
with theoretical' predictions. 

As Figures 26 through 31 show close to nadir in both polarizations, almost 
all the return comes from the large facets. But at large incidence angles- in HH 
polarization case, the- predicted return is determined by the return from medium 
size facets while in W polarization it is mainly determined by the return from small 
facets. To facilitate comparison of theoretical prediction with experimental data 
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Reflection Coefficient for Medium Size Facets for Different Moisture Contents 
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Reflection Coefficient for Small Size Facets for Different Moisture Contents 



Figure 25. Reflection coefficient for small size facets. 
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1. 2 GHz 
5.9% 
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Contribution of the Large Facets to Oq in dB 
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Contribution of the Medium Facets to Oo in dB 


Frequency - 9 GHz 
Sweep Frequency = 1. 2 GHz 
Moisture Content = 5.9% 
Polarization = VV 

S[ = 2 cm 

S| = 3 cm 



Incidence Angle In Degrees 

Figure 30. Contribution of medium size facets to the return for VV polarization. 
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(Figures 32“47), which represeni an average return over a 1 .2 GHz bandwidth, 
the theoretical calculation shown in Figures 26“47 represent the average of the 
calculated return at the indicated center frequency and the returns at the 
neighboring frequencies 400 MHz on each side . 

In fitting the theory to the data curve we use the level of the curve at 
large angles in HH polarization to determine the standard deviation of the distribution 
of the facdt size, then use the level and the shape of the curve near zero angle 
to determine the standard deviations of the slopes. Once these parameters are deter- 
mined, they are used to predict the return in VV polarization. It should be noted 
that in large angles in VV polarization the return predicted by the theory is only 
a function of the frequency and the dielectric constants of the target. 

When we increase the frequency the standard deviation of the facet size must 
stay constant, but because we include smaller facets to the large category, it is 
possible that the standard deviations of the slopes also increase. 

The results of the fittings are shown in Figures 32 through 47, where 

S = standard deviation of the slope in x direction 
= standard deviation of the slope in y direction 

S = standard deviation of the distribution of the facet size 

In order to check the accuracy of the predicted variance of the slope distribution, 
we calculate the variance of the slope of the ground from the data taken at the time 
of the radar measurement. After smoothening the surface curve. Figure 48, by the 
length of the smallest facet in tbe large catagory, i.e. X , the variance of slope was 
calculated to be 0,07, and this is very close to the predicted value of 0.075 that 
has been used to calculate the data in Figures 40 through 43, 
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Figure 32. RCS of bare ground; moisture content = 5.9%, 
frequency = 9 GHz, polarization = HH. 


no 




Figure 33. RCS of bare ground; moisture content - 5.9%, 
frequency = 9 GHz, polarization =* VV. 
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RCS of Bare Ground In dB 



Figure 34. RCS of bare ground; moisture content = 5.9%, 
frequency = 14.2 GHz, polarization = HH. 


112 



RCS of Bare Ground in dB 



Figure 35. RCS of bare ground; moisture content - 5.9 %, 
frequency * 14.2 GHz, polarization = VV. 
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RCS of Bare Ground in dB 



Figure 36, RCS of bare ground; moisture content - 9.4%, 
frequency ■ 9 GHz, polarization = HH . 


14 




Figure 37. RCS of bare ground; moisture content - 9.4%, 
frequency = 9 GHz, polarization = W. 
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Figure 38. RCS of bare ground; moisture content * 9.4%, 
frequency * 14.2 GHz, polarization = HH. 
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Figure 39. RCS of bare ground; moisture content - 9.4%, 
frequency = 14.2 GHz, polarization = W. 
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Figure 40.. RCS of bare ground; moisture content - 17.8%, 
frequency -9 GHz, polarization = HH. 
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RCS of Bare Ground In dB 



Figure 41 . RCS of bare ground; moisture content “ 17.8%, 
frequency = 9 GHz, polarization = VV. 
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RCS of Bare Ground in dB 



■ Figure 42. RCS of bare ground; moisture content - 17.8%, 
frequency * 14,2 GHz, polarization = HH . 


120 




Figure 43. RCS of bare ground; moisture content - 17.8%, 
frequency ~ 14.2 GHz, polarization =VV. 
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Figure 44. RCS of bare ground; moisture content * 21 .1 %/ 
frequency — 9 GHz, polarization = HH. 
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Figure 45, RCS of bare ground; moisture content = 21 .1%, 
frequency = 9 GHz, polarization = VV. 
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RCS of Bare Ground in dB 



Figure 46. RCS of bare ground; moisture content -21.1%, 
frequency = 14.2 GHz, polarization = HH. 
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RCS of Bare Ground in dB 



Figure 47. RCS of bare ground; moisture content -21.1 %, 
frequency ■ 14.2 GHz, polarization = VV. 
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Field No. 1 

Average Moisture Content = 5.9% 
(First 2. 5 cm) 



Field No. 2 

Average Moisture Content = 9.4% 
(First 2. 5 cm) 



Field No. 3 

Average Moisture Content = 17.8% 
(First 2.5 cm) 

Standard Deviation of Slope = 0.07 



Plate 1. Photographs of bare fields 
at four different dates. 
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Field No. 4 

Average Moisture Content = 21. 1% 
(First 2.5 cm) 



5.3 Comparison of Focet Model with Katzin's Facet Madel 

As stated earlier Kafzin assumes that the target is composed of small and 
large facets and the return power is the sum of the return from each category. 
For large facet 
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m { 




Bf 


~e^r" 


(5-13) 


where 

V = cos 0 + Z sin 0 

y 

B2. = Ka Z cos 0 

X 

cu= Kb(-Z cos 0+ sin 0) 

y 

K = 2tt/X 

X = wavelength 


P(A) 

P(Z ,Z ) 

\ y/ 


A 


1 


probability distribution function of the area of the facet, 
joint probability distribution function of the slopes in 
X and y directions 

the transition area from small to large which is equal to 

.1104x2 
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then by assuming that the slope has a Gaussian distribution; 
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277 -S* 






(5-14) 
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where 



(5-16) 


In Eq. (5-16) we can define and g 2 such that: 


-fr ^ ^ 




(5-17) 
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where 


_ tix-n frA, Z. 

<1,^- ^ 

2S)i 5^ 

^ hJt 

(j 3/7 /d^&- 

then, from Eq . (5-15) 

Az 

^ -1 T 
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(5-18) 

(5-19) 


(5-20) 


The return from the small facet is assumed to be in the Rayleigh region, hence 




•pcrtl d'fl' 


(5-21) 


where 




A 



(5-22) 


Ai and P(A) are the same as defined In Eq. (5-13). 
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Then the return from the target will be 


r = 


^5 

H 




\ rt" r>tA ) da 


(5-23) 


Now Katzin assumes that the distribution of the facet area is 


- 

A/o 


(5-24) 


If we substitute this relation in Eq . (5-23) 
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where 


Ki A 


(5-26) 


m Is the same as in Eq. (5-24). 


For the case that n 0 &n?^4,i.e. m^2 & m 4 , Eq . 5-25 becomes 
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When n = 0, then 



and when n = 4 


(5-28) 
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(5-29) 



The problem with these equations Is that they sometimes give results which are not 
acceptable. For example In Eq. 5“29 when the smallest facet of the target goes to 
zero a goes to Infinity, or In some cases one gets a negative number for cr from 
Eqs. { 5 - 27 ) and {5“28). 

One other problem is that Katzin does not take into consideration the case 
that the target is not perfectly conducting, so in order to compare his results with 
the measured data we have to only look at the shape of the curves not their level. 

In what follows we have calculated Eqs. (5~27), (5“28), and (5~29) for 
different values of the parameters; frequency, the size of largest facet, the size 
of smallest facet, shape of the facet distribution function, and standard deviations 
of the slope in x and y direction . If we compare these results as shown in Figures 49 
and 50 with the value of the measured RCS which have been shown in Figures 32 
through 47, we see that Katzin's model predicts ihe return from the bare ground very 
poorly. One important point is that for incidence angle range 20°“50° the measured 
data has a positive curvature, and this is the same as the result predicted by our model, 
while Katzin's model gives us a negative curvature; so fitting Katzin's model to the data 
is not possible. 
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RCS in dB Based on Katzin’s Derivation 
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RCS in dB Based on Katzin’s Derivation 



Incidence Angle in Degrees 


Figure 50. RCS based on Katzin's derivation. 
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6.0 CONCLUSION 


Theo^e^ical work was conducted to investigate the radar return terrain. The 
target area was modeled as a collection of facets described in terms of facet size 
and slope distribution. 

The facet model was used to calculate the radar cross section of bare ground 
and the effect of frequency averaging on the reduction of the variance of the return . 

Starting with the results obtained by Katzin [8] it was shown that when he 
divides the facets into two categories, small ancf large, he actually groups not only 
the medium size facets but also some of the small facets as large. So we divided the 
facets into three categories, small , medium and large. Katzin in his derivation 
assumes some distribution function for the area, because getting an experimental 
function for area is very difficult we changed the distribution function of the area to 
the distribution function of the length and width. We also introduced an effective 
reflection coefficient in our equations in order to be able to compare the result of the 
theory with the measured data. 

We showed that by assuming that the distribution of the slope is Gaussian, 
and by assuming that the distribution of the length of the facet is in the form of the 
positive side of a Gaussian distribution, the results are in better agreement with 
experimental data than Katzin‘s model. 

We also showed that for this calculation we do not need to know the exact 
correlation length of the small structure on the ground and introduced an effective 
correlation length which is only a function of the wavelength, i.e. the radar wave 
has .some smoothing effect. The comparison of the theory and the measured data 
shows that our conclusion is correct. 

We showed that at near zero incidence angles in both polarizations the major 
part of the return is due to large facets but at larger incidence angles it depends on the 
polarization, i.e. for HH polarization medium size facets have a major role in the 
return while in W polarization, the small facets contribute the most. 

So by knowing the level of the return at large incidence angles for HH, and 
VV polarization, and by knowing the shape and leyel of the return as a function of 
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angle near nadir we are able fo, respecHvely, calculaie shandard devia^^on of fhe 
facef size, moisfure confenf, standard deviation of the slope in y direction, and 
standard deviation in the x direction. 

In the second section of the work, namely the prediction of the number of 
independent samples due to frequency averaging, we started from Waite's, [26] 
uniform scatterer mode! and extended it and took into consideration the penetration 
effect. Comparison with data taken from alfalfa shows good agreement. We then 
assumed that the target consists of a collection of two dimensional facets and 
calculated the number of independent samples . Comparison with the data taken 
from the bare ground reveals that our result is correct for large incidence angles but 
at small incidence angles, both due to the limited number of available data points 
and the assumption of the equality of the mean and the variance in a single frequency, 
the facet model cannot predict the same number of independent samples calculated 
based on the measured data. 

We showed that based on the facet model assumption the reduction in the 
variance of the returned signal is not only a function of the product of the sweep band 
and the time span of the target, as small scatterer model indicates, but it is also a 
function of the surface properties^distributions of the slopes and facet s!ze-“, 
center frequency of the incident wave, and the polarization. 
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appendix a 


KATZIN'S ASSUMPTION 


From physical optics, the RCS for a large circular disk with a diameter D 

is given 


V- 


T- r 


2 4^ 2^{^l7/D 

xno ^ 


J 


(A-1) 


where is the grazing angle. ^ , 

Figure A-1 shows fhaf* for ycy S ^ envelope 

defined by the equation ^ • Hence its value on the average over x will 

be equal to half of the value of the envelope, i ,e, ^ j -ny? . 

So for the case when Irru^^ <9^ ^ *«®* 

(i) Is small 
(ii) p>X 

We can approximate ^ ] by ^ obtain the average RCS . 
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By substituting the relation for the area of a circular disk. 
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appendix b 


RMS HEIGHT FOR SMALL FACETS 


We would like fo calculate the r.m.s. height of small facets. Because 

from the assumption of the small perturbation technique 4 I / we assume 

that the angle is uniformly distributed between _ ^ and IL , i.e. 

H 
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-~2l JZ 


^ Fig. B-1 


(B“l) 




then because the distribution of the length of the facet is assumed to be Gaussian and 
because the length of the facets is small compared to the standard deviation of the 
facet distribution then we assume that it is uniform, i.e. 
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We know that 
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If we assume thaf 
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Figure B-3 


<io 


J p[ 

dtL 






Figure B“4 
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For ihe caselhat- U. > o ihe result wili be in the same form and the only change 
will be to substitute o( in place of -o( . So the probability density for 
will be 
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Because what- we are doing is actually integrating over shaded area in Figure B“5 
so we can change the order of integration and will have 
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If we substitute 
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appendix c 


CORRELATION RADIUS OF SMALL FACETS 


We would like to calculate the correlation radius for the small facets. 

The auto-correlation function of the surface is defined to be 

3cl+) (C-1) 



Figure C“1 

From Fig. C“1 if the distance of Z. is in a way that 't' and t-rZ 
are not on the same facet then from the basic assumption of the facet model, 
indicating that the facets are independent, we know that 

^ ® (C- 2) 
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So in order i"o contribufe f'o fhe average of Eq. (C“l) 


(C-3) 


where 


t, 


length of the facet 


a 


the distance between point t and the first point of the facet 


Because the distribution of ^ over the facet -t, is uniform, i.e. 
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then 


(C-4) 


(C-5) 


is the length of the facet and has a Gaussian distribution, i.e. 
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So by using Eqs. (C-5) and (C“6) in Eq. (C-1) 
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(C-7) 


where is the correlation d stance and ^ is the size of the largest facet in the 
small facet category. By assumption stated in Chapter 3 



So from Eq. (C~7) by putting 
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Because 



then from Eq . (C-8) 
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where 
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fhen because both A and B are small 
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From Eqs. (C-13), (C-14), (C-15) we determine 'T in a way that 
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and from Eqs. (C“13) and (C“14) 
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If we puf 



(C-20) 


then the left hand side of Eq. (C“I9) will be 

LHs ^ ^-^1 - c 


(C-21) 


so 





(C-22) 


21 In Eq. (C-22) is the effective correlation radius of the surface when it is 
looked at by a wavelength equal to X • 
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APPENDIX D 


EVALUATION OF THE INTEGRAL FOR LARGE FACETS 


In discussing the reiurn from large facefs we have to calculaie fhe integrals 
in fhe form of 
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If we intografe bofh sides vs. p fwice fhe resulf will be equal fo 
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If we expand fhe cos t'erm 
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Because fhe infegrant’ of the second integral is an odd function of u the result 
is zero. So 
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From the table [1] 
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By subsfit-uting fhe exponenHal for fhe cos funcfion 
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Eq. (D-10) becomes 
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From here 
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From riie table [1] 
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then 
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where K-| and K 2 are independent from P^or from Eq. (D-10) 
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where 



Now we will use the boundary condition principle to calculate K, and K . 

I O 
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i) ZtJO i-'O j 1 ,e . looking of a facef wUh a 

leng^h zero verHcally. In ^his case fhe relurn should be zero, 
li) K.'s.o , i .e. looking af a facet wilh 

a length zero, in an incidence angle 6. In this case the return should 
also be zero . 


From the above boundary conditions 
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and if we write this expression in the compact form 
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and from Eq. (D“I) 
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where 




(D-20) 


Now we wan^ fo calcula^e differenf'aJ'e bofh sizes of Eq. (D~l) 

vs a 
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The nexf problem is fo calculot’e u/3 • From €q, (D“3) if C is zero, fhe infegranf 
IS an odd func^ion of u so the result will be zero. But if C is not zero because 
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is very large the exponential will be very small before u gets large enough to 
have. a significant value for ^4^ pC^irVy) 

So the result of the integral will be very small, i.e . 


h/3 <>• 
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